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We investigate the statistical properties of the photon-subtractions from the two-mode squeezed 
vacuum state and its decoherence in a thermal environment. It is found that the state can be considered 
as a squeezed two-variable Hermite polynomial excitation vacuum and the normalization of this state is 
the Jacobi polynomial of the squeezing parameter The compact expression for Wigner function (WF) is 
also derived analytically by using the Weyl ordered operators' invariance under similar transformations. 
Especially, the nonclassicality is discussed in terms of the negativity of WF. The effect of decoherence 
on this state is then discussed by deriving the anal)n:ical time evolution results of WF. It is shown that 
the WF is always positive for any squeezing parameter and photon-subtraction number if the decay 
time exceeds an upper bound (kI > | In f^^)- 
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I. INTRODUCTION 

Entanglement is an important resource for quantum information [1]. In a quantum optics laboratory, Gaussian 
states, being characteristic of Gaussian Wigner functions, have been generated but there is some limitation in using 
them for various tasks of quantum information procession [Q] • For example, in the first demonstration of continuous 
variables quantum teleportation (two-mode squeezed vacuum state as a quantum channel), the squeezing is low, thus 
the entanglement of the quantum channel is such low that the average fidelity of quantum teleportation is just more 
(8 ± 2) % than the classical limits. In order to increase the quantum entanglement there have been suggestions and 
realizations to engineering the quantum state by subtracting or adding photons from/to a Gaussian field which are 
plausible ways to conditionally manipulate nonclassical state of optical field [3, In fact, such methods 

allowed the prep aration and analysis of several states with negative Wigner functions, including one- and two-photon 
Fock states fllOl flU [130, delocalized single photons lfl4l. [Tsll . and photon-subtracted squeezed vacuum states 
(PSSV). ver y si milar to quantum superpositions of coherent states with small amplitudes (a Schrodinger kitten state 
11161 [T^.[T8l. 1 19^1) for single-mode case. 

Recently the two-mode PSSVs (TPSSVs) have been paid enough attention by both experimentalists and theoreti- 
cians (T^inillllllllllllllllslillll,!!!,^ Olivares et al [11] considered the photon subtraction 
using on-off photodetectors and showed the improvement of quantum teleportation depending on various parame- 
ters involved. Then they further studied the nonlocality of photon subtraction state in the presence of noise [22]. 
Kitagawa et al [23], on the other hand, investigated the degree of entanglement for the TPSSV by using an on-off 
photondetector. Using operation with single photon counts, Ourjoumtsev et al . ifiol. ITT[| have demonstrated experi- 
mentally that the entanglement between Gaussian entangled states, can be increased by subtracting only one photon 
from two-mode squeezed vacuum states. The resulted state is a complex quantum state with a negative two-mode 
Wigner function. However, so far as we know, there is no report about the nonclassicality and decoherence of TPSSV 
for arbitrary number PSSV in literature before. 

In this paper, we will explore theoretically the statistical properties and decoherence of arbitrary number TPSSV. 
This paper is arranged as follows: in Sect. II we introduce the TPSSV, denoted as a™fe"S'2(A) |00) , where S'2(A) is 
two-mode squeezing operator with A being squeezing parameter and m, n are the subtracted photon number from 
S'2(A) 1 00) for mode a and b, respectively. It is found that it is just a squeezed two-variable Hermite polynomial 
excitation on the vacuum state, and then the normalization factor for a™&"S'2(A) |00) is derived, which turns out 
to be a Jacobi polynomial, a remarkable result. In Sec. Ill, the quantum statistical properties of the TPSSV, such 
as distribution of photon number, squeezing properties, cross-correlation function and antibunching, are calculated 
analytically and then be discussed in details. Especially, in Sec. IV, the explicit analytical expression of Wigner function 
(WF) of the TPSSV is derived by using the Weyl ordered operators' invariance under similar transformations, which 
is related to the two-variable Gaussian-Hermite polynomials, and then its nonclassicality is discussed in terms of the 
negativity of WF which implies the highly nonclassical properties of quantum states. Sec. V is devoted to studying the 
effect of the decoherence on the TPSSV in a thermal environment. The analj^ical expressions for the time-evolution 
of the state and its WF are derived, and the loss of nonclassicality is discussed in reference of the negativity of WF 
due to decoherence. We find that the WF for TPSSV has no chance to present negative value for all parameters A and 
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m, n if the decay time Kt > i In 2|±| (see Eg. ([46]) below), where n denotes the average thermal photon number in 
the environment with dissipative coefficient n. 

II. TWO-MODE PHOTON-SUBTRACTED SQUEEZED VACUUM STATES 
A. TPSSV as the squeezed two-variable Hermite polynomial excitation state 

The definition of the two-mode squeezed vacuum state is given by 

S2iX) |00) = sechAexp(a^6UanhA) |00) , (1) 

where S'2(A) = cxp[A(a^6^ — ab)] is the two-mode squeezing operator lIsoL [3ll. with A being a real squeezing 
parameter, and a, b are the Bose annihilation operators, [a, a^] — [b, 6^] = 1. Theoretically, the TPSSV can be obtained 
by repeatedly operating the photon annihilation operator a and b on 5*2 (A) |00), defined as 

|A, TO, n) =a"'6" 52(A) |00), (2) 

where |A, m, n) is an un-normalization state. Noticing the transform relations, 

Sl {X)aS2 ( A) a cosh X + b^ sinh A, 

S'|(A)&52(A) = 6coshA-l-a^sinhA, (3) 

we can reform Eq. Q as 

|A,m,n) =52(A)5|(A)a'"6" 52(A) |00) 

= 52(A)(acoshA + &^ sinh A)™ (6 cosh A -f sinh A)" |00) 

m , ,1 , 

= 52(A) sinh"+™ aV ^;^-^^6t™-;„;„tn |oo) . (4) 
^ — ' / ! m — t ! 

1=0 ^ ' 

Further note that a^" |0) = and a}\n) = ,\n-l) = 7^%a^""' |0)> leading to a'at" 1 00) ^ 

(n-;)! ^^"^' 1 00) J thus Eq.Q can be re-expressed as 

min(m,n) ; 

|A,TO,n).52(A)sinh-™A J] yT(^^?f^«^"-'^^"-' |00) 

(\ 111— I / \ n—l 

iVtanhAfen [iVtanh AaM 

^^Vr^"'""' h l\{n-inm-l)\ 

= ^ ^n+m ^2(A)g;„,„ (^VtoIhA&^ iVte^IhAan |00) , (5) 

(i%/2) ^ ' 

where in the last step we have used the definition of the two variables Hermitian polynomials ll33L[34l] . i.e., 

min(m,n) , ^fc , , „i_fe „_fe 

^ ^ (-1) m\n\(i'^ '^e" 
^m.„le,ej- 2^ fc!(m-fc)!(n-fc)! ' ^ ^ 



fc=0 



From Eq.© one can see clearly that the TPSSV |A, m, n) is equivalent to a two-mode squeezed two-variable Hermite- 
excited vacuum state and exhibits the exchanging symmetry, namely, interchanging to <^ ?i is equivalent to <^ 5^ . 
It is clear that, when to = n = 0, Eq.® just reduces to the two-mode squeezed vacuum state due to iJo,o — 1; 
while for n 7^ and to = 0, noticing i/o.ji {^,v) — y", Eq.® becomes 

C^A^o.« = n!sinh^"A), see Eq.dll]) below) 
|A, 0, n) — S2{X) \n, 0) ,which is just a squeezed number state, corresponding to a pure negative binomial state llssi] . 
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B. The normalization of |A, m, n) 

To know the normalization factor A^A.m.n of | A, m, n), let us first calculate the overlap (A, m + s,n + t |A, m, n). For 
this purpose, using the first equation in Eq.([5]) one can express |A, m, n) as 

mm(m,n) . n+m ^ i x 

I ^ X ^ / , \ m\n\ smh A coth A , , 

|A,m,n) = 52(A) ^ „ ==- \n ~ l,7n ~ I) , (7) 

which leading to 

(A, m + s, n + i |A, m, n) 



m! (n + sllJ.tsinh'^+'^+^^A 



min(m,n) , M i +i2;+s \ 

E [m + s)\n\ coih A 
n(m-/)!(n-/)!(/ + s)!' ^ 



where (5^ 4 is the Kronecker delta function. Without lossing the generality, supposing m < n and comparing Eq.® 
with the standard expression of Jacobi polynomials ll36h 



k=0 



one can put Eq.® into the following form 

(X,m + s,n + t\X,m,n) = ml {n + s)lSs^t sinh2"+" A cosh^ A/^^,"""'^) (cosh 2A), (10) 

which is just related to Jacobi poljmomials. In particular, when s = t = 0, the normalization constant Nm^n.x for the 
state \X,m,n) is given by 

Nx,m.n = {X,m,n \X,m,n) = mln! sinh^" APi""'"^^) (cosh 2A), (11) 

which is important for further studying analytically the statistical properties of the TPSSV. For the case m = n, it 
becomes Legendre polynomial of the squeezing parameter A, because of Pi*''°^(a:) = P„(x), Po{x) = 1; while for n 7^ 
and m = 0, noticing that p^^'^^ (x) = 1 then A^A.o.n = n\ sinh^" A. Therefore, the normalized TPSSV is 

||A,m,n) EE U!n!sinh2"AP^"-"'°)(cosh2A)j"^^\™6"S'2(A)|00). (12) 
From Eqs. ([2]) and (IllD we can easily calculate the average photon number in TPSSV (denoting r = cosh 2 A), 

(a^a) = 7V2„_„ msl{X)a^^+'b^-a"^+'b^S2{X) |00) 

p(n-m-l,0)/ N 

= (m+l)^2±^-— (13) 



p(n-m+l,0)/ N 

(6t5) = (n+l)sinh^A \ ^ (14) 



pr"'"'"^(T) 



In a similar way we have 



p(n-m,0)/ 

(a^b^ab) = (m + 1) (n + 1) sinh^ A 7+' ' . (15) 
Thus the cross-correlation function can be obtained by lls^. Issl. l39l] 



(2)mn _ (a^'b^'ab) 



p(n-m-l,0)/ N p("-m+l,0) / N ' 



(16) 
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Actually, the cross-correlation between the two modes reflects correlation between photons in two different modes, 

which plays a key role in rendering many two-mode radiations nonclassically. In Fig.l, we plot the graph of (A) as 

(2) 

the function of A for some different (to, n) values. It is shown that (A) are always larger than unit, thus there exist 
correlations between the two modes. We emphasize that the WF has negative region for all A, and thus the TPSSV is 
nonclassical. In our following work, we pay attention to the ideal TPSSV. 



III. QUANTUM STATISTICAL PROPERTIES OF THE TPSSV 
A. Squeezing properties 
For a two-mode system, the optical quadrature phase amplitudes can be expressed as follows: 

^ Ql + Q2 n Pi + P2 nl . ,T ry^ 

where Qi = {a + a^)/V2, Pi = [a - at)/(\/2i), Q2 = [b + fe^)/\/2 and P2 ^ {b - b^)/{V2i) are coordinate- and 
momentum- operator, respectively. Their variances are (AQ)^ = (Q^) — (Q)^ and (AP)^ = (P"^) — {P)"^ ■ The phase 
amplifications satisfy the uncertainty relation of quantum mechanics AQAP > 1/2. By using Eqs. JlQD and (1111) . it is 
easy to see that (a) = (a^) = (b) = (6^) ^ and (a^) = (a^^) = (b'^) = (6t2^ = as well as (06^) = (a+fe) = 0, which 
leads to (Q) = 0, (P) = 0. Moreover, using Eq. dlOD one can see 

, -. p(ri— m,l) / N 

(a^b^) = {ab)^^f^ -i^sinh2A. (18) 

From Eqs.dTS]), (O and it then follows that 

(AQ)2 = i((a^a) + (b^b) + {ab) + (a^b^) + 1) 

= ,p(»-L,0), M + 1) P!n+r''"\r) + in+l) PL-"'+'''\r) sinh^ A 

+ {n+l) Pi"-"'^HT) sinh2A + ^1"""'°^^)], (19) 



and 



(AP)2 = i((a^a) + (b^b) - {ab) - (a^b^) + 1) 



(„_L,o), [im + l)Pl:-r'^'\r) + (n+ l)p(r-'"+i^")(r)sinh^ A 
(n+ l)Pi"-™'i'(T)sinh2A + P/„"-™^o)(T)]. (20) 
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Next, let us analyze some special cases. When m = n = 0, corresponding to the two-mode squeezed state, Eqs. ([T9l ) 
and (I20D becomes, respectively, to 

(Ag)2| =-e^^,(^P?\ =ie-2\AQAP=i (21) 

\ \m=n=0 2 lm=n=0 2 2 

which is just the standard squeezing case; while ior m — 0, n — 1, Eqs. gi9D and fl2QD reduce to 

(AQ)^L=o..=i = (A^)l.=o,.=i = AQAP = 1, (22) 

from which one can see that the state \ \X,m,n) is squeezed at the "p-direction" when e^^'^ < i, i.e., A > iln2. In 
addition, when m — n — 1, in a similar way, one can get 

2 1 _ ^ „2\ ^ ' 2e^'*' 



lm=ri=l 2 \ 

(^^) L^«^i = 2 ^ + 2 < 2' ^23) 

which indicates that, for any squeezing parameter A, there always exist squeezing effect for state | |A, 1, 1) at the 
"p-direction". 

In order to see clearly the fluctuations of (AP)^ with other parameters m, n values, the figures are ploted in Fig.2. 
From Fig.2(a) one can see that the fluctuations of (AP)^ are always less than i when m = n, say, the state | |A, m, m) 
is always squeezed at the "p-direction"; for given m values, there exist the squeezing effect only when the squeezing 
parameter exceeds a certain threshold value that increases with the increasement of n (see Fig.2(b)). 



B. Distribution of photon number 

In order to obtain the photon number distribution of the TPSSV, we begin with evaluating the overlap between 
two-mode number state {na,nb\ and |A,m,n). Using Eq.© and the un-normalized coherent state iIstL Issll . \z) = 

exp (za^) |0), leading to {n\ — {z\\^,_^, it is easy to see that 

{na,nb \X,m,n) 

= sechA {na,nb\ a"6"e'^'^''* '^''^^ |00) 



(m + ng)! 



sechA tanh'"+"" A(5„,+,, . (24) 



It is easy to follow that the photon number distribution of | |A, m, n), i.e., 

P{na, rib) = N^X.n I ('^a, rib | A, m, n) |^ 

_ [(m -I- ria) !sechAtanh™+"° A(5„i+„^,„+„J ^ 
nJublmW. sinh^" APi"""'"^ (cosh 2A) 



(25) 
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FIG. 3: (Color online) Photon number distribution P{na,nb) in the Fock space {na,nb) for some given m — n values: (a) 
m = n = 0, A = 1, (b) m = n = 1, A = 0.5, (c) m = n = 1, A = l,(d) m = 2, n = 5, A = 1. 



From Eq. (l25D one can see that there exists a constrained condition, m+Ua — n+rib, for the photon number distribution 
(see Fig. 3). In particular, v^rhen m = n = 0, Eq. (l25D becomes 

N f sech^A tanh^"" A, Ua — rib 
U, ^^a 7^ nb 

v^rhich is just the photon number distribution (PND) of two-mode squeezed vacuum state. 

In Fig. 3, we plot the distribution P{na, m) in the Fock space ina,nb) for some given m, n values and squeezing 
parameter A. From Fig. 3 it is found that the PND is constrained by to + rta = n + nb, resulting from the paired- 
present of photons in two-mode squeezed state. By subtracting photons, we have been able to move the peak from 
zero photons to nonzero photons (see Fig. 3 (a) and (c)). The position of peak depends on how many photons are 
annihilated and how much the state is squeezed initially. In addition, for example, the PND mainly shifts to the bigger 
number states and becomes more "flat" and "wide" with the increasing parameter A (see Fig.3 (b) and (c)). 



C. Antibunching ejfect of the TPSSV 



Next we will discuss the antibunching for the TPSSV. The criterion for the existence of antibunching in two-mode 
radiation is given by ll40ll 

_(atVM^ .27) 



(at^a^) = (m + 1) (m + 2) „, , (28) 



In a similar way to Eq. ([13]) we have 



0.0 0.2 0.4 ^ 0.6 0.8 1.0 ' 0.0 0.2 0.4 ^ 0.6 0.8 1.0 

FIG. 4: (Color online) The Rab as a function of A and m, n. (1) and (2) in Fig.4 (a) denotem — and m = 1 respectively, and 
n — 2,3, 12 from down to up. 

and 

(b^^b^) = (. + !)(. + 2) sinh^ (29) 

Prn (t) 

For the state | \\,m,7i), substituting Eqs. JlSD . fl28D and ( |29D into Eq.g27D, we can recast Rat to 

jm + 1) (m + 2) Pir+T^'^°V) + (» + !)(» + 2) sinh^ APir-"+^'"' (r) ^ 
2 (m + 1) (n + 1) sinh^ AFi"-™'°)(T) 

In particular, when m — n — (corresponding to two-mode squeezed vacuum state), Eq.( |3QD reduces to Rab,m=n=o = 
— sech2A < 0, which indicates that there always exist antibunching effect for two-mode squeezed vacuum state. In 
addition, when m = n, the TPSSV is always antibunching. However, for any parameter values m, n{m ^ n), the case 
is not true. The Rab as a function of A and m, n is plotted in Fig. 4. It is easy to see that, for a given m the TPSSV 
presents the antibunching effect when the squeezing parameter A exceeds to a certain threshold value. For instance, 
when TO = and n = 2 then Rat = 6fi+yc°osh^2A) csch^ A may be less than zero with A > 0.549 about. 



IV. WIGNER FUNCTION OF THE TPSSV 

The Wigner function f WF) 1I32L I4TI, l42h is a powerful tool to investigate the nonclassicality of optical fields. Its partial 
negativity implies the highly nonclassical properties of quantum states and is often used to describe the decoherence 
of quantum states, e.g., the excited coherent state in both photon-loss and thermal channels [43^, ,44.1 . the single- 
photon subtracted squeezed vacuum (SPSSV) state in both amplitude decay and phase damping channels and so 
on llnLIillilLlii,!!!!. In tWs section, we derive the analj^ical expression of WF for the TPSSV. For this purpose, we 
first recall that the Weyl ordered form of single-mode Wigner operator [48, 49,.i50i1 , 

Ai (a) = i ■ S{a~a)5 {a* - a^) [ , (31) 

where a ^ {qi+ ipi ) / \/2 and the symbol ] ] denotes Weyl ordering. The merit of Weyl ordering lies in the Weyl 
ordered operators' invariance under similar transformations proved in Ref. ll48[] . which means 

S'' {000) [S-^ ^[S {000) S-^[, (32) 

as if the "fence" ; ; did not exist, so S can pass through it. 
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Following this invariance and Eq.® we have 

Sl (A) Ai (a) A2 (/3) ^2 (A) 

= (A) ■ 5{a-a)S{a* - a^) 6 {(3 - b) S {(3* - b^) ' S2 (A) 

= - 5 {a — a cosh A — fe^ sinh A) 5 {a* — a) cosh A — sinh A) 

(5 (/? — 5 cosh A — sinh A) 5 — 6^ cosh A — a sinh A) 

= i ■ (5 (a - a) 5 (a* - a^) (5 - 6) (5 [p* - b^) \ 
= Ai (a) A2 (^) , 

where a — a cosh A — /3* sinh X, (3 = P cosh A — a* sinh A, and p — {q2 + ip2) /\/2. Thus employing the squeezed two- 
variable Hermite-excited vacuum state of the TPSSV in Eq.® and the coherent state representation of single-mode 
Wigner operator lIsTll . 

A,{a) = e^l"!^ / ^ \z,) {-z,\e-'^^^"'-"^'^\ (33) 



where \zi) — exp (zia^ — zla) |0) is Glauber coherent state 11371 ISSfl . we finally can obtain the explicit expression of 
WF for the TPSSV (see Appendix A), 



[m!n!]^ (-tanhA)'+'-' 
l\kl [{m - ly. {n - ky.f 



\H„,-u,^k{B,A)\\ (34) 



where we have set A = — 2iQ!VtanhA, B — — 2i/3\/tanh A. Obviously, the WF W{a, (3) in Eq.(|34l) is a real function and 
is non-Gaussian in phase space due to the presence of H,n-i,n-k (B, A), as expected. 

In particular, when m = n = 0, Eq.(l34l) reduces to W{a,i3) = J^g"^'"'''"^!'^! = 
j^g2(a*/3*-(-a/3)smh2A-2(Q!Q*+/3^*)cosh2A corresponding to the WF of two-mode squeezed vacuum state; while for 

the case of m = and n ^ 0, noticing Ho,n (x, y) = and iVA,o,n = n\ sinh^" A, Eq.( |34t becomes 

W{a^P)^t^e-^\-\'~A-prL^[m'). (35) 

where i„ is m-order Laguerre polynomials and Eq.( |35D is just the WF of the negative binomial state 52(A) |n, 0) 1^35*]. 
In Figs. 5-7, the phase space Wigner distributions are depicted for several different parameter values m, n, and A. As 
an evidence of nonclassicality of the state, squeezing in one of quadratures is clear in the plots. In addition, there 
are some negative region of the WF in the phase space which is another indicator of the nonclassicality of the state. 
For the case of m and n ~ 1, it is easily seen from g35D that at the center of the phase space (a = /3 = 0), 
the WF is always negative in phase space. Fig.5 shows that the negative region becomes more and visible as the 
increasement of photon number subtracted to(== n), which may imply the nonclassicality of the state can be enhanced 
due to the augment of photon-subtraction number. For a given value m and several different values n (^ m), the WF 
distributions are presented in Fig. 7, from which it is interseting to notice that there are around \m — n\ wave valleys 
and |to — n| -f 1 wave peaks. 



V. DECOHERENCE OF TPSSV IN THERMAL ENVIRONMENTS 



In this section, we next consider how this state evolves at the presence of thermal environment. 
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FIG. 5: (Color online) The Wigner function W(a, f3) in phase space (0, 0,pi,p2) for several different parameter values m = n with 
A — 0.5. (a) m=n=0; (b) m=n=l and (c) m=n=5. 




FIG. 6: (Color online) The Wigner function W(a, /?) in three different phase spaces for m = 0, n = 1 with A — 0.3 (first row) and 
A = 0.5 (second low). 

A. Model 

When the TPSSV evolves in the thermal channel, the evolution of the density matrix can be described by the 
following master equation in the interaction picture [5Z1 

-^p(t) = (Li + L2)pW, (36) 
at 

where 

Lip = K {n + 1) {2aipa\ - a\aip - pa\a^ 

+ tin {2a\pai - aia\p - pa^a)^ , (ai = a, a2 = 6) , (37) 

and K represents the dissipative coefficient and n denotes the average thermal photon number of the environment. 
When n = 0, Eg . (1361) reduces to the master equation (ME) describing the photon-loss channel. The two thermal 
modes are assumed to have the same average energy and coupled with the channel in the same strength and have the 
same average thermal photon number fi. This assumption is reasonable as the two-mode of squeezed state are in the 
same frequency and temperature of the environment is normally the same llssl. [54ll . By introducing two entangled 




FIG. 7: (Color online) The Wigner function W(q, /3) in phase space (0, 0,pi,p2) for several parameter values m, n with A = 0.5. 
(a) ml,=n=2; (b) m=l,n=3 and (c) m=l,n=5. 
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state representations and using the technique of integration within an ordered product (IWOP) of operators, we can 
obtain the infinite operator-sum expression of density matrix in Eg .( 1361 ) (see Appendix B): 

oo 

p{t)= J2 M,,,,r,sPoMl^,^,, (38) 

i,j,r,s—0 

where po denotes the density matrix at initial time, Mij^r,s and mJ ^ ^ ^ are Hermite conjugated operators (Kraus 
operator) with each other, 



nJ + 1 y r\s\i\]\ 
and we have set T = 1 — e^^*^*, as well as 

Ti = , T2 = , T3 = ^ ' . (40) 

It is not difficult to prove the Mij^r,s obeys the normalization condition J2Tj r s=o ^1 j r s^i,j,r,s = 1 by using the 
IWOP technique. 



B. Evolution of Wigner function 



By using the thermal field dynamics theory ||55|.[56|] and thermal entangled state representation, the time evolution 
of Wigner function at time t to be given by the convolution of the Wigner function at initial time and those of two 
single-mode thermal state (see Appendix C), i.e., 

t^(a,/3,t) = ^^^-i^y ^-jpWiCv,0)e-' ™ . (41) 

Eq. (I41D is just the evolution formula of Wigner function of two-mode quantum state in thermal channel. Thus the WF 
at any time can be obtained by performing the integration when the initial WF is known. 

In a similar way to deriving Eq.( |34D . substituting Eq. (|34D into Eq. J41D and using the generating function of two- 
variable Hermite polynomials (A2), we finally obtain 



7r22"+™ {2n+lfT^D' 



-e " 



-2A-2Kt J 



EE 



f (-|tanhA)'+'= 



tt-'o nki[im-ky.in-m 



(42) 



where we have set 



C 



e 



-2Kt 



{2nT + lfC^, F 



D ' ■ / - ' ' 

/^pKt _ _ 

G = {B + B*C) , B = i2\/tanh A cosh A -I- a sinh A) , 



Ce 



Kt 



K= -j^ {A + A*C) , A = i2 Vtanh A (a* cosh X + f3 sinh A) . (43) 

Eq. (l42D is just the analj^ical expression of WF for the TPSSV in thermal channel. It is obvious that the WF loss its 
Gaussian property due to the presence of two-variable Hermite polynomials. 
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In particular, at the initial time (i = 0), noting E ^ 1, (2n + 1)^ T^D ^ 1, -I ^ 1 and ^ ^ 1, -g ^ as well as 
K A*,G^ B*, Eq.lgl) just dose reduce to Eq.(l34]), i.e., the WF of the TPSSV. On the other hand, when nt oo, 
noticing that C -> 0, £» ^ 1, S ^ 1, F ^ -1, and G/y/E 0, K/y/E -> 0, as well as iJ™,„ (0, 0) = (-1)" m\S,n,n, as 
well as the definition of Jacobi polynomials in Eq.®, then Eq.( l42D becomes 

W{a,P, oo) = -e-^^d^l'+l/^l'), (44) 

which is independent of photon-subtraction number m and ?i and corresponds to the product of two thermal states 
with mean thermal photon number n. This implies that the two-mode system reduces to two-mode thermal state after 
a long time interaction with the environment. Eq.( l44D denotes a Gaussian distribution. Thus the thermal noise causes 
the absence of the partial negative of the WF if the decay time Kt exceeds a threshold value. In addition, for the case 
of TO = n = 0, corresponding to the case of two-mode squeezed vacuum, Eq.(|42D just becomes 

VF™.„=o [a, P, t) = m-'e"*i\-\'+\mn+i^^^+"'f''\ (45) 

where 01 = tt^ (2n + if T'^D is the normalization factor, D (2n + 1)^ T^D, £ = 2 (2n + 1) T + e^^^t cosh2A, and 
2e-2''* sinh2A. Eq.lgSl) is just the result in Eq.(14) of Ref |54]. 
In Fig.8, the WFs of the TPSSV for (to = 0, rt = 1) are depicted in phase space with A — 0.3 and n — 1 for several 
different nt. It is easy to see that the negative region of WF gradually disappears as the time nt increases. Actually, 
from Eq. 043 D one can see that D > and i? > 0, so when F <Q leading to the following condition: 

1, 2f3 + 2 

Kt > Ktc= - In — — -, (46) 

2 2n -|- 1 

we know that the WF of TPSSV has no chance to be negative in the whole phase space when nt exceeds a threshold 
value ntc- Here we should point out that the effective threshold value of the decay time corresponding to the transition 
of the WF from partial negative to fully positive definite is dependent of to and n. When nt ~ ntc, it then follows from 
Eq.(l42l) that 

W(a,P,U = tanh"+"Asech^A ^_,2...[|^|2^|g|2_(^,^,^^g^t^„i^^] 



,;a (^i Vtanh A/3*e''*'= , iVtanh Aa* e*^* 



(47) 



which is an Hermite-Gaussian function and positive definite, as expected. 

In Figs. 9 and 10, we have presented the time-evolution of WF in phase space for different n and A, respectively. 
One can see clearly that the partial negativity of WF decreases gradually as n (or A) increases for a given time. This 
case is true for a given n (or nf) as the increasement of nt (or n) . The squeezing effect in one of quadratures is shown 
in Fig. 10. In principle, by using the explicit expression of WF in Eq.( l42D . we can draw its distributions in phase space. 
For the case of to = 0, rt = 2, there are two negative regions of WF, which is different from the case of to = 0, n = 1 
(see Fig. 11). The absolute value of the negative minimum of the WF decreases as nt increases, which leads to the full 
absence of partial negative region. 



VI. CONCLUSIONS 



In summary, we have investigated the statistical properties of two-mode photon-subtracted squeezed vacuum state 
(TPSSV) and its decoherence in thermal channelwith average thermal photon number n and dissipative coefficient 
K. For arbitrary number TPSSV, we have for the first time calculated the normalization factor, which turns out to 
be a Jacobi polynomial of the squeezing parameter A, a remarkable result. We also show that the TPSSV can be 
treated as a squeezed two-variable Hermite polynomial excitation vacuum. Based on Jacobi polynomials' behavior 
the statistical properties of the field, such as photon number distribution, squeezing properties, cross-correlation 
function and antibunching, are also derived analytically. Especially, the nonclassicality of TPSSV is discussed in terms 
of the negativity of WF after deriving the explicit expression of WF. Then the decoherence of TPSSV in thermal 
channel is also demonstrated according to the compact expression for the WF. The threshold value of the decay time 
corresponding to the transition of the WF from partial negative to completely positive is presented. It is found that 
the WF has no chance to present negative value for all parameters A and any photon-subtraction number (m, n) if 
Kt>^ In for TPSSV. The technique of integration within an ordered product of operators brings convenience in 
our derivation. 
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FIG. 8: (Color online) The time evolution of WF {m = 0,n — 1) at (gi,q2,0,0) phase space for n = 1, A = 0.3. (a) nt = 0.05, (b) 
Kt = 0.1, (c) Kt = 0.12,(d) Kt = 0.2. 




FIG. 9: (Color online) The time evolution of WF (m = 0,n — 1) in (qi, g2, 0, 0) phase space for A = 0.3 and Kt = 0.05 with (a) 
n = 0,(b) n = l,(c) n = 2,(d) n = 7. 




FIG. 10: (Color online) The time evolution of WF (m = 0,n = 1) in (gi,g2,0,0) phase space for n — l,and Kt — 0.05 with (a) 
A = 0.03,(b) A = 0.5,(c) A = 0.8,(d) A = 1.2. 
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FIG. 11: (Color online) The time evolution of WF for m = 0,n = 2 in (gi,g2,0,0) phase space with (a) Kt — 0,(b) K,t — 0.05, (c) 
Kt = 0.1, (d) Kt = 0.2. 
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Appendix A: Deriviation of Wigner function Eq.([34l) of TPSSV 

The definite of the WF of two-mode quantum state j^I') is given by W{a, (3) — Ai (a) A2 (/3) |^), thus by uisng 
Eqs.®, and ([33]) the WF of TPSSV can be calculated as 

W{a, (3) = (A, m, n\ | Ai (a) A2 (/?) | |A, m, n) 
sinh"+'" 2A 



, ^ , , {00\H,n.n (-iVtanhAfe, -iVtanh \aj Ai (a) 

A2 (/3) Hr,^^n (iVtanh Afe^ iVtanh Aa^) |00) 



sinh"+"' 2A 



kll -N2r-2(ziQ*-azJ)-2(z2/3*-/3^5) 



X i7„ 



^— i\/ tanh Az2, — tanh Azi^ Hm,n ^— iV tanh Azj , — i V tanh AzJ" 
Further noticing the generating function of two variables Hermitian polynomials, 



exp [-tt' + ei + rt'] 



Eq. (AI) can be further rewritten as 



sinh"+"2A ^2|a|^+2|,3|' 5"+" 5"+" 



-tT-t'r' 



-|zi|^ + (-2a*-iVtanh AT)zi + (2a- j\/tanh At')z* 



7 ^2 



t=T=0 



(Pz2 _ 



|z2|^ + (-2/9*-iVtanh At)z2 + (2/3-iv'tanh At')z* 



sinh"+"'2A _„^,2_,isP 9™+" 9™+" 



-2|a| = -2|/3|^ 



X e 



tT-t'T'+A*T/+B't'+AT+Bt-(u'+TT') tanh A 



t=T=t'=r'=0 



where we have set 



(AI) 



(A2) 



(A3) 



B = -2i/3VtanhA, A = -2ia:VtanhA, 



(A4) 



14 



and have used the following integration formula 



-ie-^,Re(C) <0. 



Expanding the exponential term exp [— {tf + tt') tanh A] , and using Eq.(A2), we have 

sinh"+"2A ^ ^ ..-.u ^^l+k 



W{a,(3) 



Ql+k Ql+k 



«,2 (-tanhA) 
;=0 k=o 



g-2|a|^-2|^|' 

Q2m 



Ukl 



Q2n 



-tT+AT+Bt-t'r'+A'Ti+B't' 



sinh"+" 2A 



t=T=t'=r'=0 
l+k 



-2|c5|=-2|^|' (-tanh A) 

Ukl 



1=0 k=0 



Ql+k 



Ql+k 



:Hm,n {B,A)Hm,n {B* , A*) . 



dB^dA^ dB*^dA*k 
Noticing the well-known differential relations of Hm,n (e, e) , 

, ^ m\n\H,n-i,n-k{€,e) 
ae'aefc ""'"^"'^^ {m-l)\{n-k)\ ' 

we can further recast Eq.(A6) to Eq.( |34D . 

Appendix B: Derivation of solution of Eq. (l36b 

To solve the ME in Eg . (1366 . we first introduce two entangled state representations Js 



1 77a) = exp 
\T]b) = exp 



-\\vb\'' +vbb^ -r]iy^ +b''y^ 



1 00) 
00), 



which satisfy the following eigenvector equations, for instance. 



{a~d)) \'qa) = \ria) , {a'^ - a) \r]a} = vl \Va} , 

{l]a\ (a'f - a) = 77* {rial , {Va\ (« " O^) = Va {Va\ ■ 



(A5) 



(A6) 



(A7) 



(Bl) 
(B2) 

(B3) 



which imply operators (a — d'^) and (a^ — a) can be replaced by number rja and 77*, [(a — a'^), (a'^ — a)] = 0. Operating 
two-side of Eq. (l36D on the vector \Ia,Ib) = Iva = 0) ® |7yb = 0), (denote \p{t)) = p{t) /(,)), and noticing the 
corresponding relation: 



a\Ia,Ib) = a'^ \Ia,Ib) , \IaJb) = d\Ia,Ib) , 
b\Ia,Ib)^P\IaJb), b^\IaJb)^b\Ia,Ib), 



(B4) 



we can put Eq. ([36]) into the following form: 



— \p {t)) = [k {n + 1) (2aa — a — a^a) + nn (2a^a^ — aa^ — aa^) 



+ K{ 



; (n -I- 1) (2bb - bh - Pbj + Kfi (2b^P - bb^ - W 

It's formal solution is given by 

\p (t)) = exp ^Kt {n + 1) (2ad — a^a — a^a) -|- Ktn (2a^d^ — aa^ — aa^) 
+Kt in + 1) (2bb - b^b - b^b^ + ntn (2b^b^ - bb^ - bb^^ \po) , 



(B5) 



(B6) 
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where |po) = po \Ia,Ib)- In order to solve Eq.(B6), noticing that, for example, 

2aa — a) a — a) a = — (a) — a) (a — a^) + aa — a) a\ 

we have 

\p (t)) = exp [(aa — a) a) + l) Kt\ 



(B7) 



X exp 
X exp 
X exp 



2n+ 1 



(1 - e^^*) (at - a) (a - S^) 



(^66 - 6^6^ + 



2n+ 1 



(fot - &) (fo - 6^) 



Ipo) 



where we have used the identity operator, exp[A(A + aB)] — e^^ exp[cri?(l — e ^'^)/t] valid for [A, B] = tB. 
Thus the element of p (t) between (ria,r]b\ and \Ia,Ib) is 



{Va,Vb\ p{t)) ^ exp 



— ^1 ('7a| 



('7a 



(B8) 



(B9) 



from which one can see clearly the attenuation due to the presence of environment. 

Further, using the completeness relation of \ria, rjb), J Iva, Vb) {Va,rib\ = 1 and the IWOP technique ll58l,[59ll . 



we see 



(fria(frib 



Va,Vb) {Va,Vb\ P{t)) 



{nT+l) 



X exp 
X exp 



a^a + + + Vb^ InTa 



T3 \ aa + bb 



Pa \Ia,Ib) , 



(BIO) 



where Ti,T2 and T3 are defined in Eq.( l40t . Noticing Eq.(B4), we can reform Eq.(BlO) as p{t) — 

J2Tj,r,s=o ^^i.j,r,sPoMjj^^ ,,,where Mij,r,s and M/j_r,s are defined in Eq.([39l). 
Appendix C: Deriviation of Eq. fl41D by using thermo field dynamics and entangled state representation 
In this appendix, we shall derive the evolution formula of WF, i.e., the relation between the any time WF and the 

initial time WF. According to the definition of WF of density operator p: W (a) = Tr [A (a) p], where A (a) is the 

single-mode Wigner operator, A (a) = {2a) (—1)° By using {n \m) = Sm.n we can reform W (a) as [i6Q1 



W{a) = V(n,n|A(a)p|m,m) = - (C=2a| p) , 



(CI) 



where (^| is the conjugate state of whose overlap is (77! ^) = i exp [i (^77* — Cv)] ^a Fourier transformation kernel. 
In a similar way, thus for two-mode quantum system, the WF is given by 



W {a, P) = Tr [A, («) Afc (/3) p] = — (ea=2a, 6=2/3! p) ■ 
Employing the above overlap relation, Eq.(C2) can be recast into the following form: 

Then substituting Eq.(B9) into Eq.(C3) and using the completeness of (^|, / ^ |^) (^| — 1, we have 



X e 



e ''\vbe-^'\ 



Po) 



(fVa(fVb 2^T\{^^\^ + \nb?) 



X e 



a'ria—avl+P'Vb — PVb 



{Vat 



47^2 

, Vbe~ 



6=2C;6=2r;) 



(C2) 



(C3) 



(C4) 
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Performing the integration in Eq.(C4) over (PrjadPrn, then we can obtain Eq.( l41D . 

Making variables replacement, > C,, * Eg. ( 141 D can be reformed as 

W {a, f3, t) = 4e4«* J d^Cd^^W^f (C) Wt iv) 

xTy|e''*(a~\/rc),e"*(/3-\/T7;),o}, (C5) 

where W*^ (C) is the Wigner function of thermal state with average thermal photon number n: W^'^ {Q — 

TT(2n+i) ^ 2*^+1. Eq.(C5) is another expression of the evolution of WF and is actually agreement with that in 
Refs.ll53l,[54ll. 
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